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ABSTRACT 

The  temperature  distribution  due  to  a  rotating  or  dithering  Gaussian  laser  beam  on  a  finite  body  is  obtained 
numerically.  The  authors  apply  various  techniques  to  solve  the  nonhomogeneous  heat  equation  in  different 
spatial  dimensions.  The  authors  ’approach  includes  the  Crank-Nicolson  method,  the  Fast  Fourier  Transform 
(FFT)  method  and  the  commercial  software  COMSOL.  It  is  found  that  the  maximum  temperature  rise  de¬ 
creases  as  the  frequency  of  the  rotating  or  dithering  laser  beam  increases  and  the  temperature  rise  induced 
by  a  rotating  beam  is  smaller  than  the  one  induced  by  a  dithering  beam.  The  authors  ’numerical  results  also 
provide  the  asymptotic  behavior  of  the  maximum  temperature  rise  as  a  function  of  the  frequency  of  a  rotating 
or  dithering  laser  beam. 

Keywords:  Crank-Nicolson,  Fast,  Fourier  Transform  (FFT),  Nonhomogeneous  Heat  Equation,  Rotating 

or  Dithering  Gaussian  Laser  Beam 


INTRODUCTION 

The  impact  of  laser  beams  on  metals  or  other 
materials  is  of  great  interest  in  industry  and 
military.  This  is  mainly  due  to  the  fast  processing 
time  and  precise  operation.  Adeeper  understand¬ 
ing  of  the  physics  requires  a  modeling  of  the 
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heating  process.  There  are  many  previous  works 
on  the  theoretical  and  numerical  modeling  of 
temperature  profdes  induced  by  laser  radation 
in  solids  (Araya  &  Gutierrez,  2006;  Bertolotti 
&  Sibilia.,  1981;  Burgener  &  Reedy,  1982; 
Calder  &  Sue,  1982;  Cline  &  Anthony,  1977; 
Lax,  1977  and  1978;  Moody  &  Hendel,  1982; 
Sanders,  1984).  However,  mostofthese  studies 
are  limited  to  a  scanning  Gaussian  beam.  Until 
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Figure  1.  (1)  A  dithering  laser  beam  on  a  1-D  rod.  (2)  1-D  temperature  distribution  along  the 
rod  from  various  numerical  methods.  (3)  ID  maximum  temperature  rise  of  steel  AISI 4340  ver¬ 
sus  frequency  of  the  dithering  laser  beam.  (4)  The  curve  in  (3)  is  well  approximated  by  the 
function  =  137.0864/frequency  +  51.6960  ■ 


Dithering  Laser  Beam  g(x,i)=  /„e 


(1) 

600 


(3) 


recently  laser  forming  of  plates  using  a  rotat¬ 
ing  or  dithering  laser  beam  has  been  studied 
(Sistaninia  et  ah,  2009).  A  more  detailed  study 
of  the  temperature  rise  induced  by  a  rotating  or 
dithering  laser  beam  on  a  semi-infinite  domain 
is  newly  provided  by  Zhou  (20 1 1).  In  this  paper 
we  want  to  extend  the  study  in  Zhou  (2011) 
to  a  more  realistic  finite  geometry  and  figure 
out  the  quantitative  relationship  between  the 
maximum  temperature  rise  and  the  frequency 
of  the  rotating  or  dithering  laser  beam. 

It  should  be  pointed  out  that  the  models 
used  here  are  heat  equations,  which  are  deter- 


1-D  Methods  Comparison 
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ministic.  In  order  to  be  able  to  include  random 
effects  (e.g.  laser  beam  jitter)  it  is  desirable  to 
use  stochastic  differential  equations  instead  of 
heat  equations.  This  will  be  the  direction  of 
future  work. 

We  organize  our  paper  into  five  sections. 
In  the  first  four  sections  we  present  the  numeri¬ 
cal  modeling  of  the  temperature  distributions 
induced  by  a  dithering  or  rotating  laser  beam 
in  one-dimensional,  two-dimensional  and 
three-dimensional  finite  solids,  respectively. 
Conclusions  and  future  work  are  given  in  the 
last  section. 
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1-D  MATHEMATICAL 
MODELING 

Consider  a  laser  beam  hitting  a  1-D  rod  with 
finite  length  L  .  The  beam  moves  along  the  rod 
which  is  insulated  at  the  two  endpoints.  Math¬ 
ematically,  the  temperature  distribution  of  the 
rod  can  be  modeled  by  the  nonhomogeneous 
heat  equation: 


du  d\i 

dt  ~  dx"^ 


(1) 


The  boundary  conditions  reflect  the  as¬ 
sumption  that  no  energy  escapes  into  the  ambient 
at  the  air/rod  interface.  This  is  a  good  approxi¬ 
mation  for  most  materials  under  consideration 
because  heat  flow  by  conduction  through  the 
material  is  much  larger  than  heat  loss  by  radia¬ 
tion  or  convection  at  the  air/material  interface. 

An  analytical  solution  based  on  the  ei¬ 
genfunction  expansion  can  be  derived  for  this 
initial-boundary  value  problem.  Assume 

(4) 


where  u  (x,  denotes  the  temperature  rise  of 
the  rod  at  position  x  (o  <  x  <  ij  and  time  t , 
ttj.  is  the  thermal  diffusivity  of  the  rod,  is 

the  thermal  conductivity,  and  q[x,t^  is  the 
energy  distribution  of  the  moving  laser  beam. 
In  the  case  of  a  dithering  laser  beam  shown  in 
Figure  1  (1),  q[x,t^  can  be  expressed  as 


I  \  ■  2-11? 


(2) 


1  tr  f  \  f  1  r,  \ 

where  X^yxj  =  cos  -  \n  =  1, 2, •  •  ■  j 


are 


the  eigenfunctions.  After  substituting  (4)  into 
(I),  one  has 


n2 


(5) 


with  the  initial  condition  (oj  =  0  .  Equation 

(5)  is  a  first-order  linear  ODE  for  and  it  can 
be  solved  by  the  standard  integrating  factor 
method: 


where  x^  is  the  position  of  the  dithering 
Gaussian  beam,  x^  is  the  initial  position  of  the 
laserbeam,  is  the  intensity  ofthe  laser  beam, 
is  the  effective  radius  of  the  laser  beam,  and 
A:  is  a  constant  used  for  the  Gaussian  model. 
The  initial  condition  for  m(x,  tj  is  zero  which 
assumes  that  the  rod  has  the  same  temperature 
as  the  ambient  initially.  The  boundary  conditions 
impose  that  the  rod  is  insulated  at  the  two  ends : 
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Once  the  coefficients  (t  j  s  are  found, 
the  solution  u[x,t^  is  given  analytically  by  (4). 
Flowever,  the  computation  of  (4)  can  be  expen¬ 
sive  if  one  uses  direct  summation.  A  more  ef¬ 
ficient  numerical  approach  is  to  use  the  Fast 
Fourier  Transform  (FFT).  In  order  to  so,  one 
first  needs  to  even  extend  the  problem  from  the 
domain  [O;^]  to  and  make  it  periodic 

with  period  2L  .  Then  one  can  apply  FFT  and 
its  inverse  to  find  the  nuiuerical  solution.  Details 
of  the  implementation  can  be  found  in  Tan’s 
thesis  (Tan,  20 lO). 
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An  alternative  way  to  solve  the  non-ho- 
mogeneous  1-D  heat  equation  (1)  with  insulat¬ 
ing  boundary  conditions  is  to  apply  the  Crank- 
Nicolson  method.  We  begin  by  discretizing  the 
interval  [o,  i]  into  n  subintervals  with  nodes 


1 

i - 

2 


Ax, 


A  L 

Ax  =  — , 


i  =  1, 2,  •  •  • , « 


(7) 


Then  we  approximate  (1)  with  the  Crank- 
Nicolson  scheme: 


jt+l  k  I  k  f\  k  t  k  i+l  fA  i+l 

If,  —  U,  CL^  I  —  2u,  +  U,_,  ^  V 


—  2w.*  +  u 


At 


Ax- 


Ax^ 


+  +  '  =  1,2, ■■■,«;  4=iAr 

(8) 


and  the  insulating  boundary  conditions  can  be 
approximated  by 


u  ^  =  u 

n+1  n  ’ 


for  all  k  >  1. 


(9) 


In  (8)  M.*  =  u{^iAx,kAt^ .  The  linear 
system  (8)  and  (9)  can  be  solved  directly  using 
MATLAB  built-in  function  or  the  Thomas  al¬ 
gorithm  (Burden  &  Faires,  2005). 

A  third  way  to  solve  (1)  is  to  use  the  com¬ 
mercial  software  COMSOL.  COMSOL  Multi¬ 
physics  is  a  finite  element  analysis  software  for 
various  physics  and  engineering  applications. 

In  Figure  1(2)  we  plot  the  numerical  solu¬ 
tion  of  (1)  obtained  from  these  three  different 
methods  at  t  =  .  As  shown  in  Figure  1(2), 

the  temperature  rise  along  the  rod  from  differ¬ 
ent  methods  agrees  with  each  other  very  well. 
The  two  sharp  peaks  of  the  temperature  rise 
occur  near  x  =  0.25  and  x  =  0.75 ,  which  are 
the  endpoints  of  the  path  of  the  dithering  laser 
beam  given  in  Equation  (2).  Since  the  beam 
moves  with  the  prescribed  path 
x^  j  =  X|j  +  a  sin  2TTty/ T  where 

Xp  =  0.5,  a  =  0.25,  T  =  1 ,  the  velocity  of  the 
beam  is  x^  (t)  =  0.5Ttcos(2TTtj .  Around  the 
two  endpoints  sin^2TTtj  =  ±l  and 


Table  1.  Material  properties  of  Steel  AlSl  4340  and  parameters  of  the  1-D  dithering  laser  beam 


Property 

Name 

Value 

Unit 

Heat  Capacity 

C 

P 

475 

J/(kg*K) 

Density 

P 

7850 

Kg/(m^3) 

Thermal  Conductivity 

Kr 

44.5 

W/(m*K) 

Thermal  Diffiisivity 

“r 

Kt 

m''2/s 

Melting  Point 

1783 

K 

Magnitude  of  Gaussian  beam 

1.0e9 

W/(m''3) 

Effective  radius  of  Gaussian  beam 

0.02 

m 

Dithering  frequency 

1 

f 

1 

Hz 

Center  of  dithering 

0.5 

m 

Dithering  radius 

a 

0.25 

m 
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cos(2tt^)  =  0  .  Consequently,  the  beam  moves 
slower  and  stays  longer  near  these  two  end¬ 
points.  Thus,  the  temperature  is  higher  near  the 
two  endpoints.  The  material  that  we  choose 
here  is  steel  AISI  4340  with  its  physical  prop¬ 
erties  listed  in  Table  1  and  the  parameters  of 
the  dithering  laser  beam  are  also  given  in  Table 
1. 

Figure  1(3)  depicts  the  maximum  tem¬ 
perature  rise  as  a  function  of  the  frequency  of 
the  dithering  beam  at  a  fixed  time  t  —  0.1 5  . 
Figure  1(3)  shows  that  the  maximum  tempera¬ 
ture  rise  acts  as  a  decreasing  function  of  the 
frequency  of  the  dithering  beam.  In  other  words, 
the  faster  the  dithering  beam  moves,  the  lower 
the  maximum  temperature  rise.  When  the  beam 
does  not  dither,  the  maximum  temperature  rise 
can  be  as  high  as  534if  .  Once  the  beam  dithers, 
the  maximum  temperature  decreases  quickly. 

The  asymptotic  behavior  of  the  maximum 
temperature  rise  can  be  predicted  numerically 
by  a  least-squares  curve  fitting.  In  Figure  1(4) 
we  show  that  the  curve  in  Figure  1(3)  can  be 
well  approximated  by  a  function 

^max  “  137. 0864y/ frequency  4-  51.6960  ,  Our 
result  indicates  that  at  a  fixed  time  the  maximum 
temperature  rise  decays  asymptotically  as  a 
reciprocal  function  of  the  frequency  of  the 
dithering  laser  beam  as  the  frequency  goes  to 
infinity.  An  intuitive  physical  explanation  goes 
like  this.  After  long  time  the  temperature  rise 
attains  pseudo-steady  state.  If  we  observe  the 
temperature  rise  at  a  fixedpoint,  the  temperature 
rise  oscillates  periodically  as  the  beam  passes 
through  this  point  periodically.  So  the  period 
of  the  temperature  rise  oscillation  coincides 
with  the  period  of  the  beam.  The  average  tem¬ 
perature  rise  over  one  oscillation  period  is  in¬ 
dependent  of  the  frequency  of  the  beam.  The 
maximum  temperature  rise  over  one  oscillation 
period  minus  the  average  temperature  rise  is 
proportional  to  the  period,  which  is  inversely 
proportional  to  the  frequency  of  the  dithering 
laser  beam. 


Now  we  conduct  an  asymptotic  analysis 
for  an  infinite  long  rod  and  show  that  the  lead¬ 
ing  order  term  of  the  maximum  temperature 
rise  indeed  is  proportional  to  1/frequency  as 
the  frequency  of  the  dithering  beam  gets  large. 
For  simplicity,  consider  the  following  1-D  non- 
homogeneous  heat  equation 


dll 


a » ^  -^0 

dx 


2'Tr/| 

X  —  a  sin 

—  oo  <  :\:  <  CX3 


(10) 


with  initial  condition 

It  (x,  oj  =  0  (— cxD  <  X  <  ooj .  The  solution  of 
(10)  can  be  found  by  Duhamel’s  principle 
(Asmar,  2004;  John,  1981): 


-n. 


—  exp 


dyds. 


(11) 


We  need  to  borrow  a  well-known  result 
from  probability  theory  (Hayter,  2002)  to 
simplify  (11). 

Lemma.  Suppose  and  are  nonzero 
constants.  Then 


f  — 

-oo 


;exp 


(  \2 

/  \2l 

[x-y] 

1 

(y-^) 

Ac: 

1 - 

Ac: 

cfy 


pexp 


+< 


(12) 


In  other  words,  the  integral  ofthe  product  of 
two  Gaussian  functions  turns  out  to  be  another 
Gaussian  function.  The  lemma  expresses  the 
fact  that  the  convolution  of  two  Gaussian  (or 
normal)  densities  is  another  Gaussian.  The  proof 
of  the  lemma  is  straightforward  and  is  skipped 
here.  A  simple  2-D  version  of  the  proof  can  be 
found  in  Zhou  (2011). 
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Applying  the  lemma,  we  arrive  at 


oo 

/ 


exp 


4(t-5) 


exp 


y  —  asm - 

T 


4Tt  7  -  s  +  ft; 


■exp 


4c 


.  2t\s 

X  —  a  sin - 

T 


dy 


Ait -s  +  r‘ 


(13) 


Substituting  (13)  into  (11)  yields 


0  yjt-s  +  r^ 


exp 


.  2'K5 

X  —  a  sin - 

T 


4(7 -.  +  r;) 


ds. 


(14) 


Makingachangeofvariables  =  t  —  s 
and  then  dropping  the  subscript,  we  obtain 


u{xj]  =  21, 


exp 


X  —  a  sin  z 


ds 


(15) 


J  [oj  cos  (fe)  -I-  sin  (fc)], 

k=0  ^ 


(17) 


where  and  are  the  Fourier  coefficients 
that  depend  on  x,  s  and  a  . 

Using  (17),  (15)  becomes 


00 

k=l 


(18) 


where 


II 

ds, 

cos  k 

+  b.  sin 


/t  >  0. 

(19) 


Note  that  the  leading  term  is  which  is 
O  (1)  order,  independent  of  e  .  We  continue  to 
estimate  U,  . 

k 

When  Q ,  there  is  no  singularity  in  the 
integrand  function  in  .  Integration  by  parts 
yields 

U,^0(z).  (20) 


where 


211(7  —  s] 

(.)  =  4(.  +  r;), 

e  =  T,  z=  ''  ' 

e 

(16) 


We  wish  to  find  the  asymptotic  expansion 
of  It  (x,  7 j  for  small  values  of  e  .  To  do  so, 
consider  the  Fourier  series  expansion  of 

—  flsinzj 


as  a  function  of  z  : 


So  asymptotically  m(x,  tj  behaves  as 
O  (1  j  +  O  (e  j  when  ,  which  is  exactly 

what  we  have  observed  in  F igure  4  even  though 
the  rod  considered  there  is  finite  instead  of 
infinite. 

It  is  worthwhile  to  point  out  that  if  =  0  , 
then  there  is  a  singularity  in  the  integrand  func¬ 
tion  in  .  Applying  the  standard  method  of 
stationary  phase  (Bush,  1992)  as  in  Zhou  (20 11), 
one  can  show  that  ui^x,t^  then  behaves  as 

o(l)  +  o(x/IJ  as  e  — >  0  ,  which  shares  the 

same  scaling  law  as  the  3-D  semi -infinite  case 
(Zhou,  2011). 
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Figure  4.  (1)  Snapshots  of  the  normalized  temperature  rise  for  steel  AISI  4340  induced  by  a 
rotating  Gaussian  beam  where  a.^.  =  1.19  x  10“^.  (2)  2D  maximum  temperature  rise  versus  fre¬ 
quency  of  the  rotating  laser  beam  for  a  thin  steel  AISI  4340  film.  (3)  The  numerical  results  in 
(2)  fitted  by  function  =  293.5153/frequency +  50.8908- 
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2-D  MATHEMATICAL 
MODELING 

When  a  laser  beam  shines  on  a  two-dimension¬ 
al  rectangular-shaped  thin  film,  the  temperature 
rise  of  the  film  u  (x,  y,  is  governed  by 


du 

-  =  OL^ 

dt  ^ 


dx^ 


d"u 

d/ 


0  <  X  <  L  <  y  <  L  , 


(21) 


tion  expansion  method  and  end  up  with  double 
summation,  which  is  computationally  expen¬ 
sive.  An  alternative  approach  is  to  use  Fast 
F ourier  Transform  (FFT) .  Again,  one  first  needs 
to  even  extend  the  problem  from  the  rectangu¬ 
lar  domain  |o,  L  to  a  larger  domain 

\—L  .L  \  x  \  —L  and  then  apply  FFT  and 

L  X'  X  \  I  y  \  r  r  j 

its  inverse  to  obtain  the  numerical  solution. 

Another  straightforward  method  to  solve 
(21)-(24)  is  to  apply  the  Crank-Nicolson 
scheme: 


where  the  heat  source  /(x,  due  to  a  rotat¬ 
ing  laser  beam  can  be  written  as 


/  \  2Trt 

XcV)  =  x„  +  acos  — 


T 

y,(t)  =  y„+.sin^, 


4 

L 


(22) 


In  the  case  that  b  —  0  in(22)thelaserbeam 
is  called  a  dithering  beam.  The  2-D  rotating 
laser  beam  and  the  finite  thin  film  are  illus¬ 
trated  in  Figure  2(1)  while  the  heat  source 
f{x,y,  t  j  given  in  (22)  is  shown  in  Figure  2(2). 

The  boundary  conditions  impose  the  insu¬ 
lating  condition: 


dll  „  T 

—  =  0  at  X  =  0  or  L  , 
dx 

—  =  0aty=0ori  , 
dy 


(23) 


—  2u.  ,*  +  u.  j 

‘  -  2u  +  11.  , 

2 

—  2u.  +  u.  , 

/  -  2«. 

1  1  i.j+i  i,j  i,]-i 

2 

A>.‘ 

*  A/ 

(25) 


where  the  grid  points  are  given  by 


1 

i - 

2 


Ax, 


L 

Ax  = 
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1 

'“2 


Ay,  Ay  =  ^, 
m 


yj  = 

=  kAt,  k  =  0,1,2, 


i  =  1, 2,  •  •  • ,  H 
j  =  l,2,---,m 


(26) 


and  M. =u\x., y. , j .  The  boundary  condi¬ 
tions  are  approximated  by 


=  u. 


n+lj 

k 

,m+l 


(27) 


and  the  initial  condition  is 

w(^,>’,0)  =  0.  (24) 

Like  1-D  case,  one  can  derive  the  analyti¬ 
cal  solution  to  the  initial-boundary-value 
problem  (IB  VP)  (2 1  )-(24)  using  the  eigenfunc- 


The  linear  system  (25)  together  with  (27) 
can  be  solved  efficiently  in  MATLAB  by  tak¬ 
ing  advantage  of  the  built-in  sparse  matrix 
operations. 

A  third  method  to  solve  the  2D  problem 
(21)-(24)  is  to  apply  COMSOL  directly. 

Before  we  present  our  numerical  results, 
we  would  like  to  investigate  the  effects  of 
physical  parameters  on  the  solution  of  (21).  In 
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Figure  2.  (1)  A  schematic  2-D  diagram  shows  a  rotating  laser  beam  and  the  finite  work  place. 
(2)  The  plot  of  the  heat  source  due  to  a  Gaussian  beam. 
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order  to  do  so,  we  first  rescale  the  time  variable 
as  ~  ^  and  then  drop  the  subscript.  Then 

(21)  can  be  rewritten  as 


dll 

d\i  d\ 

1  ^ 

t 

-  = 

-  — 1”  - 

“1 - f 

x,y, — 

dt 

dx^  dy^ 

k/ 

“rj 

(28) 


Substituting  (22)  into  (28),  we  obtain 


dll 

dt 


d\i  du 
dx^ ""  d/ 


Let  u  denote  the  solution  to  the  following 
IB  VP  problem 


dt 


^1 


+  exp 


dx 

du 

dy 


0  at  X  = 
0  aty  = 


0  or  L^, 
0  orL^. 


Then  it  is  obvious  that  the  solution  to  (28) 
with  same  initial/boundary  conditions  is  related 
to  the  solution  of  (30)  through  the  relationship: 

u  (x,  y,  t)  =  ii  (x,  y,t\.  (31) 

ZTta  Aj, 

So  from  a  computational  point  of  view, 
one  should  solve  (30)  to  obtain  the  normalized 
solution  u{x,y,t^  and  then  multiply  the  nor- 

malizedsolutionbyaconstantfactor 
to  get  the  solution  for  (28),  or  equivalently  (21). 
From  this  simple  analysis,  it  is  clear  that  the 
effect  of  the  thermal  diffusivity  a^,  is  to  change 
the  effective  frequency  of  rotation  and  the  effect 


of  the  intensity  of  the  laser  beam  and  the 
thermal  conductivity  is  to  rescale  the  tem¬ 
perature  distribution  by  a  factor  proportional 
to  ^0  / ^7-  .  The  same  argument  applies  to  the 
ID  and  3D  cases. 

Now  we  show  some  numerical  results.  In 
Figure  3  we  take  the  snapshots  of  the  numerical 
solutions  of  the  problem  (30)  with  a^,  =  0.1 
at  various  times  where  the  parameters  for  the 
rotating  laser  beam  are  listed  in  Table  2.  Our 
numerical  experiments  have  shown  that  the 
solutions  obtained  from  the  Crank-Nicolson 
method,  FFT  and  COMSOL  all  agree  with  each 
other  (Tan,  2010),  which  we  do  not  present 
individually  here.  We  only  give  the  FFT  results 
in  Figure  3.  One  can  see  from  Figure  3  that  the 
maximum  temperature  rise  increases  as  time 
increases.  The  hottest  spot  is  the  place  that  is 
directly  hit  by  the  laser  beam. 

Figure  4(1)  shows  the  numerical  results  of 
(30)  forsteelAISI4340 where  a^.  =  1.19  x  10“® 
(refer  to  Table  1).  Due  to  the  small  valueof  a^, 
heat  does  not  diffuse  fast  enough  and  conse¬ 
quently  a  hot  ring  starts  to  form.  Using  the 
values  listed  in  Table  2,  we  find  that  the  normal¬ 
ized  maximum  temperature  in  Figure  4(1  )d 
corresponds  to  a  value  270.6794Ai  after  mul¬ 
tiplying  by  the  factor  ■  Similarly, 

the  corresponding  time  of  Figure  4(1  )d  is 
0.0977^  after  multiplying  by  . 

In  Figure  4(2)  we  plot  the  maximum  tem¬ 
perature  rise  for  different  values  of  the  fre¬ 
quency  (i.e.  the  reciprocal  of  the  period  T )  at 
a  fixed  time  t  =  0.1^  for  a  thin  steel  film  with 
the  shape  of  a  unit  square.  As  the  frequency 
increases,  the  maximum  temperature  rise  de¬ 
creases  and  the  asymptotic  behavior  of  this 
decrease  is  depicted  in  Figure  4(3)  by  a  fitting 
function 

=  293.5153/frequency  +  50.8908  .  Fig¬ 
ure  4(3)  implies  that  the  maximum  temperature 
rise  is  proportional  to  the  reciprocal  of  the 
frequency  as  the  frequency  gets  large. 
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Figure  3.  Snapshots  of  the  temperature  rise  induced  by  a  rotating  Gaussian  beam  at  various 
times  where  ctj.  =  0.1 


In  Figure  5(1)  we  give  the  snapshots  of  the 
temperature  rise  when  the  laser  beam  is  dither¬ 
ing  (where  a  —  0.25  and  b  —  0  '\n  equation 
(22)).  Here  a^,  =  0.1  and  the  solutions  are 
obtained  by  solving  (30)  with  FFT.  Note  that 
the  dithering  beam  moves  at  its  slowest  speed 


at  the  endpoints  of  its  trajectory  so  these  two 
endpoints  are  exposed  directly  to  the  beam 
longer  than  other  places. 

In  Figure  5(2)  we  compare  the  maximum 
temperature  rise  against  the  frequency  of  the 
rotating  and  dithering  beams  for  a  square-shaped 
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Table  2.  Parameters  of  the  2-D  rotating  laser  beam 


Parameter 

Value 

Unit 

“r 

1 

ni^ls 

Kr 

1 

W/(m*k) 

T 

l.OeS 

Wlm^ 

d 

0.02 

m 

1 

m 

L 

y 

1 

m 

X, 

0.5 

m 

a 

0.25 

m 

h 

0.25 

m 

period 

1 

s 

thin  steel  film  at  t  =  0.1^ .  As  the  frequency 
increases,  the  maximum  temperature  rise  de¬ 
creases.  The  maximum  temperature  rise  cor¬ 
responding  to  the  rotating  beam  drops  much 
faster  than  the  one  induced  by  the  dithering 
beam  for  the  same  frequency.  It  implies  that 
rotating  a  beam  is  more  effective  than  dithering 
a  beam  if  one  desires  a  lower  temperature  rise. 


3-D  MATHEMATICAL 
MODELING 

In  the  case  of  a  rectangular  prism  hitting  by  a 
laser  beam,  the  temperature  rise  on  the  solid 
body  is  described  by 


d'^u  d'^u  d'- 

<  X  <  L  <  y  <  L  ,0< 
—  —  ^  — 


du 

—  =  a 

dt 

0  < 


'll 

<  z  <  L_ 


q{x,y,z,t), 


(32) 


where  the  heat  source  q[x,y,z,t^  due  to  the 
rotating  laser  beam  can  be  modeled  as 

q  (^,  y,  z,t)  =  f  (x,  y,  t'j  6  (z) .  (33) 


Here  /(x,  y,  t  j  is  given  in  equation  (22) 
and  the  delta  function  in  z  expresses  the  as¬ 


sumption  that  all  the  energy  is  absorbed  at  the 
surface  ^  =  0  which  is  hit  directly  by  the  beam. 
The  boundary  conditions  are  insulating  and  the 
initial  condition  is  zero. 

Our  numerical  experiments  in  ID  and  2D 
cases  have  indicated  that  COMSOL  yields  reli¬ 
able  results.  So  for  3D  case,  instead  of  using 
FFT  or  the  Crank-Nicolson  method,  we  use 
COMSOL  alone  to  compute  the  temperature 
rise  induced  by  a  rotating  laser  beam.  The 
parameters  of  the  3D  rotating  laser  beam  are 
specified  in  Table  3.  Details  on  how  to  imple¬ 
ment  COMSOL  can  be  found  in  Tan  (2010). 

Figure  6(1)  gives  several  snapshots  of  the 
temperature  rise  induced  by  a  rotating  Gaussian 
laser  beam  at  different  times  within  one  time 
period.  The  hottest  spot  is  the  place  that  is  hit 
directly  by  the  laser  beam,  as  in  the  2D  case.  The 
maximum  temperature  rise  reaches  1468A^  at 
time  t  Those  points  far  away  from  the 
heat  source  have  little  temperature  rise.  Figure 
6(2)  is  similar  to  Figure  6(1),  but  for  a  different 
period  0.1^  .  Figure  6(2)  shows  a  maximum 
temperature  rise  drops  from  ITOSAT  with  period 
1^  to  670K  with  period  0.1^  ,  which  is  about 
54%  reduction  in  the  maximum  temperature 
rise.  Therefore,  the  maximum  temperature  rise 
can  be  significantly  reduced  by  increasing  the 
frequency  of  the  rotating  laser  beam. 

In  Figure  7(1)  we  show  the  temperature 
rise  at  different  layers  at  time  t  =  1^  corre- 
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Figure  5.(1)  Snapshots  of  the  temperature  rise  induced  by  a  dithering  Gaussian  beam  at  various 
times  with  =  0.1 .  (2)  Comparison  of  the  2D  maximum  temperature  rise  versus  frequency 
of  the  rotating  and  dithering  laser  beams  for  a  square-shaped  thin  steel  fdm. 
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Table  3.  Parameters  of  the  3-D  rotating  laser  beam 


Input 

Value 

Unit 

h 

5.0e5 

Wlni^ 

d 

0.02 

m 

4 

1 

m 

L 

y 

1 

m 

\ 

0.5 

m 

A 

0.5 

m 

a 

0.25 

m 

h 

0.25 

m 

period 

1 

s 

spending  to  the  results  in  Figure  6(1).  It  is 
observed  that  heat  does  not  spread  out  downward 
quickly  and  there  is  almost  no  temperature  rise 
for  the  body  part  that  is  more  than  0.1m  below 
the  top  surface.  This  is  because  the  diffusivity 
of  steel  AISI  4340  is  quite  small  (recall 
ttj.  =  1.19  X  10“^ )  and  thus  heat  diffuses  very 
slowly.  For  materials  with  larger  diffusivity  we 
observe  that  the  heat  spreads  out  fast  and  the 
temperature  rise  occurs  in  other  layers  below 
the  top  surface. 

Finally,  the  relationship  between  the  maxi¬ 
mum  temperature  rise  and  the  frequency  of  the 
rotating  laser  beam  is  depicted  in  Figure  7(2). 
As  expected,  the  increase  of  the  frequency  of 
the  rotating  laser  beam  leads  to  a  decrease  of 
the  maximum  temperature  rise.  In  Figure  7(3) 
we  fit  the  numerical  results  in  Figure  7(2)  by 
a  function 


T 


1228.3 


frequency 


522.5 


(34) 


CONCLUSION  AND 
FUTURE  EFFORTS 

We  have  calculated  the  temperature  rise  induced 
by  a  rotating  or  dithering  Gaussian  laser  beam 
for  a  finite  body  in  different  spatial  dimensions. 
We  have  confirmed  that  for  a  finite  solid  the 
maximum  temperature  rise  can  be  reduced 
by  increasing  the  frequency  of  a  rotating  or 
dithering  beam.  Furthermore,  the  temperature 
rise  induced  by  a  rotating  laser  beam  is  smaller 
than  the  one  induced  by  a  dithering  laser  beam. 
We  have  given  the  asymptotic  behavior  of  the 
maximum  temperature  rise  as  a  function  of 
the  frequency  of  the  dithering  or  rotating  laser 
beam.  Future  efforts  will  include  (1)  experi¬ 
mentally  determining  the  actual  temperature 
profile  induced  by  a  rotating  or  dithering  laser 
beam  on  different  materials  and  comparing  the 
measured  temperature  profiles  to  the  predicted 
profiles;  (2)  using  stochastic  differential  equa¬ 
tions  to  revisit  the  problem;  (3)  investigating 
the  temperature  rise  in  a  two-layer  structure. 


using  the  lease-squares  approach.  The  good 
agreement  between  the  fitted  curve  and  the 
numerical  results  indicates  that  the  maximum 
temperature  rise  is  inversely  proportional  to  the 
frequency  of  the  rotating  beam. 
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Figure  6.  Snapshots  of  the  temperature  rise  on  a  steel  AISI  4340  box  induced  by  a  rotating 
Gaussian  laser  beam  using  COMSOL  (1)  with  period  Is;  (2)  with  period  0.1s 
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Figure  7.  (1)  Temperature  rise  of  a  box  made  of steel  AlSl  4340  at  different  depth  corresponding 
to  Figure  6(1):  (a)  ^  =  1  (top  surface  hit  directly  by  the  laser  beam)  (b)  z  =  0.99  (c)  z  =  0.95 
and  (d)  ^  =  0.90  .  (2)  3D  maximum  temperature  rise  of  steel  AISI 4340  versus  frequency  of  the 
rotating  laser  beam.  (3)  The  numerical  results  in  (2)  are  fitted  by  the  function 

=  1228.3/frequency  +  522.5  ■ 
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